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Abstract 

We consider the ultra-discrete Burgers equation. All variables of the equation are discrete. We classify 
the equation into five regions in the parameter space. We discuss behavior of solutions. Using this 
equation we construct the deterministic surface growth models respectively. Furthermore we introduce 
noise into the ultra-discrete Burgers equation. We present the automata models of the KPZ equation. 
One model corresponds to the discrete version of the ASEP and the other to the Kim-Kosterlitz model. 
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1 Introduction 



The noisy Burgers equation appears in a variety of problems in non-equilibrium statistical mechanics, jjj 
Denoting the velocity field by u(x,t) the equation reads 

du d 2 u „ du ,„ „ . 

m = d^ + 2u d- x +n01se - (L1) 

Burgers studied the noiseless equation with random initial data. We are interested in the case the noise 
is of the form d^/dx, to conserve u locally. Then ([O]) is a prototype for a driven diffusive system. If we 
set h = J udx, (Id) becomes 

dh d 2 h ,dh. , „, 

m = o^ + ^ + ^ (L2) 

We can interpret h(x,t) as the height of a (one dimensional) surface. (|l.2| ) is called the Karder-Parisi- 
Zhang (KPZ) equation. Q This equation governs the shape fluctuations of the various growth model. || 

In this paper we investigate the Burgers equation and the surface growth model without noise and 
with noise. With noise, the Burgers equation describes steady growth in the long time limit, however, 
without noise, it becomes the relaxation of an initially rough surface to the flat surface. 

The relations between the Burgers equation and the surface growth model as studied from the view 
point of the universality classes, because the surface growth models are the discrete and the Burgers 
equation is continuous. We study the direct relations between discrete and continuous using the ultra- 
discrete method. Recently Nagatani, Nishinari and Takahahi presented the ultra-discrete version of the 
Burgers equation We study solutions of the this equation. Behavior of solutions can be classified 

into five regions in the parameter space. The solutions can be understood from the viewpoint of creation 
and annihilation of particles and anti-particles. Using behavior of particles and anti-particles we can 
construct the surface models. The surface model corresponds to the deterministic version of the KPZ 
equation. Further more we present ultra-discrete version of the noisy Burgers equation and the KPZ 
equation. We show the direct relations between the cellular automata and KPZ equation. One model of 
these cellular automata is the time discrete version of the asymmetric simple exclusion process (ASEP) 
which belongs to KPZ universality class. The other model corresponds to the restricted solid on solid 
(RSOS) model which is introduced by Kim and Kosterlitz. || 

This paper is organized as follows. In the section 2 we classify the equation into five regions in 
the parameter space and discuss the behavior of the solutions from the viewpoint of particles and anti- 
particles. In the section 3 we study the symmetry of the equation. We explain why the automata model 
represent the property of continuous Burgers equation. In the section 4 we discuss the stability of the 
equation. We study the meaning of the ultra-discrete limit. In the section 5 we construct the deterministic 
surface growth models. In the section 6 we introduce the noise into the discrete Burgers equation. In 
the section 7 we construct the noisy surface growth models using the ultra-discrete KPZ equation. The 
models are equivalent to the well known models which belong to the KPZ universality. The last section 
is devoted to the concluding remarks. 

First we review the derivation of the ultra-discrete Burgers equation. Discretizing of both time and 
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space variables of diffusion equation, we can obtain 



-i 



(1.3) 



AT {AX) 2 

where AT and AX are lattice intervals in x and t respectively. Using the discrete analogue of Cole Hopf 
transformations 

h 

where c is constant. We can obtain the lattice version of the Burgers equation 



where <5 = AT/(AX) 2 . 

We introduce a parameter e and new variables 



4=e^ /e , ^—^ = e- M '*, c 2 = e L / e . (1.6) 
J CO 



Then (1.5) becomes 



E/ t+1 = [/' + elogl + expH )+cx P H 2±i )] 

e e 

- M XJ\ i+Ut -L 
- elog[l + exp(^— ) + cxp( J - 3 )] (1.7) 

Here we take the so called "ultra-discrete limit" || 

e — >0 + . (1.8) 

Then we can obtain 

U!j +1 = U) + min(M, V]_ Xl L - Uj) - mm(M, U$, L - U j+1 ). (1.9) 
Here we use the relation 



A B , 



lim elog(e« + e~) =max(A,B). (1.10) 



We call (1.9) the ultra-discrete Burgers equation. 



2 Classification 



(1.9) has a solution 



U] = - + max(0, K(j + 1) + fii + 9 ) - max(0, K(f) + nt + 6 ), (2.1) 



with 



\ max(i - M, K, -K) - § + M for M < |, 

where Oo is constant. This solution corresponds to the famous shock wave solution of the Burgers 
equation. Fig .|l| 
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Figure 1: Shock wave solution 



Here we classify five cases. 

(Al) 

(B) 
(CI) 



\ + \K\ < M, (A2) \<M<\ 

— = M, 
2 



\K\ 



- 2 >M>I-\K\ 



(C2) - - \K\ > M. 



(2.3) 



At first we study the case (A) and (B). Note that in this case the dispersion relation is = \K\. The 
sign of K chooses the propagating direction and the velocity of the wave is —1 (+1) for K > (K < 0). 



If we set XJ\ > § , then (O) becomes C/ t+1 = Uh x . On th other hand, setting Uj < (p~9|) becomes 



( "j—v I n both cases equations become linear. Then stable rectangular wave becomes a solution. 



We can write the rectangular solution explicitly as 



U 



L 



- + max(0, K(j + 1) 



Qt + 9 ) - max(0, K(j) +Slt+ 9 ) 
-max(0, A(j + 1 + 1) + Sit + 9 ) + max(0, K(j + l) + Qt+ 6 ). 



(2.4) 



The size of the rectangular is K x Z Fig. 2. 



L/2 











K 











L/2 



Figure 2: Rectangular solution 

Here we introduce an interpretation of particles and anti-particles. We call that particles (anti- 
particles) are created at time t and j when Uj > -| (Uj < -j). Uj — -j can be understood as the number 
of particles or an anti-particles. In this mean this particles or anti-particles are Boson. Then we call 
Uj = as vacuum state. From this interpretation this equation has chairarity. Particles go only to 
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the left and anti-particles go only to the right. If there are only particles (anti-particles), creation and 



annhiration of particles (anti-particles) do not occur. It can be seen from the result (2.4) is stable solution 



of 0. 

We need to study interaction of particles and anti-particles. To consider interaction we consider 
collision of particles and anti-particles. Annihilation of particles and anti-particles occurs. If same 
numbers of particles and anti-particles meet, they become vacuum. The configuration of particles and 
anti-particles are in Fig|| If the number of particles (anti-particles) are larger than the number of 
particles of anti-partials (particles), then only particles (anti-particles) whose number is difference of 
particles and anti-particles survive. 



t 



t + 1 



t + 2 



t + 3 



Figure 3: Collision of particles and anti-particles 
Note that the point where particles and anti-particles meet does not move. It can be seen also in 



continuous limit as the collision of the shock waves. In the Burgers equation ( |1 . l|) setting u(x,t) = 
u(x — vt), where v is propagation velocity. Imposing the boundary condition u —> u± and du — > for 
x — > ±oo, then we can obtain the soliton condition 

u + + u_ = — v (2-5) 

In the limit u + = it_ this condition implies v — 0. It means that the collision point does not move in 
the continuous Burgers equation. In fact the number of particles (anti-particles) is larger than that of 
anti-particles (particles), the collision point moves to the left (right) and the particles (anti-particles) 
only survive. It can be seen in the soliton condition ( |2.5| ). 

The other configuration is in Fig. 4. In the case (A) creation of pairs of particles and anti-particles 
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between particles and anti-particles occurs. On the other hand in the case (B) particles and anti-particles 
part over without creation of pairs of particles and anti-particles. It is the crucial difference between (A) 
and (B). 

We consider the case same number K of particles and anti-particles part over. In the case (Al) K 
particles and K anti-particles are created. On the other hand in the case (A2) creation of M — j pairs 
of particles and anti-particles occurs. If the number of particles are K\ and anti-particles are K 2 , in the 
case (Al) K pairs are created, in the case (A2)creation of M — -| pairs occurs, where min(ifi, K 2 ) = K. 
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(A) 



t + 1- 



t + 2- 



(B) 



t + 1- 



t + 2 



Figure 4: Separations of particles and anti-particles in the cases (A) and (B) 
Next we consider the case (C). We can rewrite the dispersion relations 

(CI) = K + M - ^, (C2) n = 0. (2.6) 



The characteristic is the dispersion relations with a gap in the case (CI). From (|2.6|) in the case (C2) the 



shock wave solution (2.1) does not move. In other words Uj 1 ' = Uj . Then any initial configurations 



of particles and anti-particles does not move. In the case (CI) the speed of the shock wave solution is 

\K\+M-% , . 

2 -, (2.7) 



K 

where v is the speed of the shock wave. The shock wave with positive parity, K > 0, propagates left 
with negative velocity, whereas the shock wave with opposite parity, K < 0, propagate in the forward 
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direction as we have shown in the cases (A) and (B). However in the case (CI) the speed v has a range. 
The maximum speed is 1 or — 1. Notice that in the limit K — > ±00 the speed of the shock wave become 
±1. If the shock waves with the positive parity and negative parity collide, the annihilation of the particles 
and anti-particles occurs as we have seen in the cases (A) and (B). The difference is the speed of the 
annihilation. In the case (CI) there are not rectangular solutions. Then the pair creation of particles and 
anti-particles does not occur. From these results the case (CI) represents the property of the continuous 
Burgers equation. 

3 Symmetry of Equations 



Here we note that the Burgers equation (1.1) is invariant under the parity transformation x — ► —x 
provided u — > —u. This feature is related to the presence on a single spatial derivative in the derivative 
term. We can not see this parity invariant in the discrete version of the Burgers equation. We consider 
the other ultra-discrete limit 

e — > -0. (3.1) 

In this limit we use a following relation 

lim elog(e^ + e"?) = mini A, B). (3.2) 

Using this relation we can take the other ultra-discrete limit. Then we can obtain 

L/j +1 = C/j + max(M, U]_ x , L -Uf) -mBx(M,U*,L -Uij +1 ) 

= XJ) + min(M, U] +1 ,L - U}) - min(M, Uj, L - U^), (3.3) 



where L = L and M = L — M. (3__3) can be obtained using the transformations of (JL 



j + j — 1 — ► j + 1, M — > M, L^L. (3.4) 



This is the parity transformation. (3.4) corresponds to x — » —x in the continuous case. (3.1) is the 



transformation Uj — ► — Uj. Then the invariance under parity transformation is recovered. The particles 
move only to the right and the anti-particles only to the left for ([T^). The particles and anti-particles 
move to the opposite directions to (fT 



Here we consider why the parity invariance can be seen in the continuous case (1.1) and ultra-discrete 



limit ( |l.9| ), but not be seen in the full discrete type euqtaion (1.5). The relation between the discrete 
variables and the continuous variables are following 

1 u' 

u(jAx,tAT) = — log^. (3.5) 
Aa c 



Using (1.8), we can obtain a relation 



u(jAx,tAT) = -^(Uj - ~). (3.6) 
The vacuum state Xl\ — L/2 corresponds to u = 0. It can be seen in the collision of the shocks. 
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Using the relation (3.E) and (3.6) the relation among continuous, discrete and ultra discrete variables 



arc 



v ~ log u\ - U] 



(3.7) 



Then the the parity invariance is the same in the continuous (1.1) and ultra-discrete Burgers equations 



1.9) 



4 Stability 



Here we consider stability of the solution of the discrete Burgers equation (|L5|). We substitute plane 



wave solution u* = e^ k+UJt into the discrete diffusion equation (F3). We can obtain the condition for the 
stability condition |e w | < 1, 

5 = 0. (4.1) 



Using the second and the third equations of ( [Dj| ) we can get the relation 

6 = 



2 + exp[i(f -M)] 

Using this relation, the cases (A), (B) and (C) correspond to 6 — 1/2, (5=1/3 and (5 = respectively.Fig.5. 



(4.2) 



From the condition (4.1) the cases (A) and (B) are unstable and (C) is only stable in the discrete Burgers 
equation. It is the reason why the case (C) represents the property of the continuous Burgers equation 
(See section 2). It also can be seen that we need to set <5 — ► to obtain the continuous Burgers equation 



(1.1) from the lattice Burgers equation (F5). On the other hand there is the rectangular solutions in the 
cases (A) and (B) which can not be seen in the Burgers equation. The rectangular solutions of (A) and (B) 
correspond to solutions of configurations of two solitons with opposite parity in the continuous Burgers 
equation. To construct this configuration we need the condition that the solitons are well separated 
and non overlapping. Then the size of the solutions are oo x oo. They are rescaled and becomes the 
rectangular solutions of the cases (A) and (B) in the ultra discrete limit. 
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Figure 5: Stability of discrete Burgers equation 



5 Deterministic ultra-discrete KPZ equation 

If we set u(x,t) — (d/dx)h(x,t), h(x,t) is governed by the deterministic KPZ equation. To construct 



automata model for deterministic KPZ equation, from (3.7) we can set 

L 



The last term is selected the constant for the vacuum state. We rewrite (1.9) 



L 



H] +1 = max^.i, B) + -- M, H* j+1 ) + C. 



where C is a constant. In the same way we can rewrite (3.2) by the independent variables H 



Hf 1 = minCfll.! , H) + - - M, H] +l ) + C, 



(5.1) 



(5.2) 



(5.3) 



where C is a constant. We call these equations (5^) and fl5.3| ) "ultra-discrete deterministic KPZ equa- 
tion." 

To describe interface dynamics we consider the periodic condition 



rrt rrt 

H 1 - H N . 



(5.4) 



Here Hf denotes the height of the surface at the lattice site i at the integer time t. Solutions of ultra- 
discrete KPZ equation are classed by (A1)-(C2) as the ultra-discrete Burgers equation. 

We consider the deterministic dynamical model. In this interface model the integer height variable 
may differ on neighbor sites j and j + 1 only ±1. The two elementary steps, deposition and evaporation, 
which define the surface evolution are illustrated in FigJ|. In a lattice gas language the height differences 
between neighorbouring sites are mapped to a particle occupation number n x = 0, 1 with the presence of 
a particle on x corresponding to slope —1 between sites j — 1 and j in the interface model and a vacancy 
at site x corresponding to slope +1. The particles have a simple dynamics. With rate 1 they jump to the 
right (left) except when the final site is occupied, in which case they stay (hard core exclusion). In this 
model growth (evaporate) can occur in local minima (maxima). These models are deterministic. The 
configuration at time t + 1 is determined from the configuration at time t by a local rule which depends 
on nearest neighbors only. The corresponding lattice gas evolves according to the automaton rule 184. Q 

In the case (Al) we set K = ±1, C = —1, C = 1, in fl5.2| ) and (5.3) describe following deterministic 
model, 

H] +1 = max( J ffj_ 1 , H] +1 ) - 1, (5.5a) 



and 



^ +1 =min(^_ 1 ,i^ +1 ) + l. 



(5.5b) 



The model particles move only to the right corresponds to ( 5.5a ) and the model only to the left is 



(5.5b). The above deterministic dynamical model is equivalent to the ultra-discrete KPZ equation (5.5a) 
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Figure 6: Deterministic Surface Growth model A 



or ( 5.5b ). ( 5.5a ) is the surface evaporate model and ( 5.51: ) is the surface growth model. We call these 
models "A model" . 

We consider another interface model. In this interface model the integer height variable ffj may differ 
on neighbor sites j and j + 1 only ±lor 0. Consider a one dimensional surface configuration parallel to 
the j axis of a square lattice. 



In the case (B) we set K = ±1,0, C = 0, in (5.2) describe following deterministic model, 

H) +1 =max(ff*_ 1 ,iTj, J ffj +1 ). 
The deterministic model (3) was studied in M. 



(5.6) 



In the case (B) we set K = ±1,0, C = -1 and C = 1, in (^2|) and (|5.3|) 

H*+ x =wax{H t j _ 1 ,H t j ,H t j+1 )-l, 



and 



H) +l =min(H t j _ 1 ,H t j ,H t j+1 ) + l. 



(5.7a) 



(5.7b) 



In the model (5.7a) a particle is evaporated at the local maxima and in the model (5.7b) a particle is 



deposited at the local minima. The dynamical equations (5.7a) and (5.7b) restrict the height difference 



±1 or automatically. We call the models (5.7a) and (5.7b) " B model ". 

In the case (C2) we also consider as the RSOS model. If we we set K = ±1, 0, C 



-1 and C = 1, 



in ( |5.2| ) and ( |5.3| ), we can obtain the model whose dynamics is trivial H* +l = Hj + 1 or = H j — 1 
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Figure 7: Deterministic Surface Growth model B 

6 Ultra-discretization with noise 

Next we consider the discrete Burgers equation with noise, 



1 ,MJ ± l u t 



a j U J 1 



1-2(5 1 ,,t , 1 t( 



(6.1) 



o5 



where £* is the noise. If we set £* = 0, (6.1) is the discrete Burgers equation (1.5). In the limit AX — > 



and AT — > with 5 — > 0, (O) becomes the noisy Burgers equation 



du d 2 u ^ ® u d£_(x,t) 
dt dx 2 dx dx 1 

where the white noise £(x,t) 

< i{x, t)£(x',t') >= BS(x, x')S(t, t'). 

The relation between the discrete noise and the continuous noise is 

1 



We introduce a parameter e and new variables 



(Axy 



We take the ultra-discrete limit ( |1.8| ) in the discrete noisy Burgers equation fl6.1| ) . We can obtain 
Uj +1 = U] + max(-M, -C/J, C/j - L, - -) - max(-M, -£/•_!, t/,- - L, X d - -). 

7 Ultra-discrete KPZ equation 



(6.2) 



(6.3) 



(6.4) 



(6.5) 



(6.6) 



We rewrite this equation using the independent variables H l - (5.1) 



H) +1 = maxCHj.i, Hj +-- M, H] +1 ,H] + X}) + C, 



(7.1) 



where C is is a constant. This is the noisy version of the model (p.2[). In the same way we can obtain 



= nrinCfll,!, H} + -- M, Hlj +1) H] + If) + (5, 



(7.2) 
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where C is is a constant. This is the noisy version of the model (5.3). 
Here we set K = ±1, C\ — — 1, and C2 — 1, we can obtain 

Hp 1 = max^, ffj + X*) - 1, 

and 

H] +1 = min(fZ|_ 1 ,fl1 +ll H) + Xj) + 1. 



(7.3a) 



(7.3b) 



These are the noisy version of the A models ( |5.5a ) and ( 5.5b| ). Here we set the noise as following 



\rt \T t 



1 probability q 
— 1 probability 1 — q 



(7.4) 



In ( 7.3a ) and ( 7.3b ) we consider the non-deterministic dynamical model. In the lattice gas language 
the particles have a simple dynamics. With rate q they jump to the right except when the final site is 



occupied for ( 7.3a ). For ( 7.3b ) with rate q they jump to the left. In ( 7.3a ) for t time steps the probability 
distribution that a particle moves n space steps is the binomial distribution 

t 



P{t,n) = 



q n (i- q y 



(7.5) 



It is well known that in the large t limit with A = qt =const, we can obtain the Poisson distribution. 
Then in the continuous limit of time we can obtain the Poisson process. 

A lattice model of particles moving stochastically with hard-core exclusion is called the asymmetric 
simple exclusion process (ASEP).|l Each particle hops to the right (left) nearest site with the probability 



Pi&t (p2<it) in every infinitesimal time interval dt. It is a Poisson process. (7.3a) corresponds to the time 



discrete version of ASEP with p\ — and (7.3b) corresponds to the ASEP with p^ = 0. 
For the B model we set K = ±1,0, C\ = —1, and C2 — 1, we can obtain 



m +i 



maxiH^^HlH^Hj + Xj) - 1, 



and 



H] +1 = min^.i, Hi, H f j+1 , Hi + Xf) + 1, 



(7.6a) 



(7.6b) 



In the model (7.6a) a particle is evaporated at the local maxima with the rate q and in the model (7.6b) 
a particle deposited at the local minima with the rate 1 — q. 

In the RSOS model introduced by Kim and Kosterlitz (KK),|| a particle is deposited at randomly 
selected site as long as the height difference Ah between nearest-neighbor columns remains as h < 0. 
Then a particle is deposited in the local minima. The difference between our B model and the KK model 
is time steps. In the B model in one time step m particles arc deposited. It corresponds to m time steps 
in KK model. The difference does not affect the universality. 



13 



8 Concluding Remarks 



We consider the ultra-discrete Burgers equation. We classify the equation into five regions in the pa- 
rameter space. We discussed the behavior of the solutions. The A and B models have the rectangular 
solutions. On the other hand the C model has only shock wave solutions. We construct the determin- 
istic surface growth models using A and B models. Furthermore we introduce the noise into the ultra 
discrete Burgers equation. We present the automata models of KPZ equation. The A model becomes the 
discrete time version of ASEP and B model corresponds to the KK model. It is well known ASEP and 
KK model belong to the universality of KPZ equation. We discussed the symmetry of the ultra-discrete 
and continuous Burgers equation. We showed why automata model represent the property of continuous 
Burgers equation. Furthermore we study the stability of equations. Using this result we can understand 
the parts of the meaning of ultra-discrete limit. 

We can extend these method to the higher dimension and we can introduce the quenched noise into 
the ultra-discrete KPZ equation. About these contents we will discuss else where. 
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